QCD analysis of xFg at NNLO: 
the theoretical uncertainties 
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Abstract 

The next-to-next-to-leading order (NNLO) QCD analysis of the experimental xF^ 
structure function from CCFR data is performed. The theoretical uncertainties of 
the analysis are discussed. 



The study of DIS structure functions has been 
a fruitful source of information on the structure 
of the proton and for testing perturbative QCD. 
Between the most outstanding results of the 
program of analyzing structure functions in QCD 
was the early explanation of the scaling violation 
phenomena. In the last two decades the second 
order approximation (NLO) has been extensively 
compared with F2 and xF^ data. More recently 
the third order terms (NNLO) have been fully 
calculated for the case of the coefficient functions m 
but only partially for the anomalous dimension part 
(n = 2,4,6,8 for singlet and non-singlet operators, 
and n = 10, only in the non-singlet case) B. This 
has allowed the analysis of xF^ y, Q and F2 (non 
singlet H and singlet g) at NNLO. 

In this note we review the theoretical uncertain- 
ties involved in the analysis of the structure function 
xFq at NNLO. For that task, we firstly remind the 
method of calculation and the most relevant results 
of the fits. 

The QCD evolution of the moments is given by: 



Af„(g2) /a,(Q2) 






AD{n,Q^)C(^\Q^) 
M^iQl) ~ \as{Ql) ) AD{n,Ql)C(^){Ql) 

where AD and C come respectively from the 
anomalous dimensions and coefficient function 
terms in the renormalization-group equation (see 
the explicit forms in M). The running coupling 
constant as{Q^) is obtained from the expression 
in terms of inverse powers of IniQ^ / K\-—). Target 
mass corrections are also added in the calculation 
of the moments (see [H). 



The moments in Eq. (1) at the initial scale are 
MniQl) = So dxx'^-'^Ax^{l - xfil + 72:). 

The structure function is reconstructed from 
its moments by using the expansion in terms of 
orthogonal Jacobi polynomials: 



N. 
n=0 



n 



xF, = x'^il-xf V e^^''(x)V4"'(a,/?)M,+2(Q') 






(2) 

where c" {a, 0) are combinatorial coefficients given 
in terms of Euler F-functions of the a and [3 weight 
parameters which have been fixed to 0.7 and 3 
respectively by the reasons given in [Q,^]. 

Power corrections are included in the analysis 
using two different approaches. Firstly, in the form 
given by the Infrared Renormalon Model (IRR) 
g adding in Eq. (1) the contribution M^^^ — 

C{n)Mn{Q'^)A2/Q^, with A2 a free scale parameter. 
Secondly, adding in the r.h.s. of Eq. (2) the term 
h{x)/Q'^ , with h{x) a free parameter for each a;-bin 
of the data set. 

Table 1 summarizes the results of the fits to XF3 
CCFR data |9|. For comparison we have applied 
the same kinematic cuts as in Ref. B, i.e. Q^ > 5 
GeV2, x < 0.7 and W^ > 10 GeV^. At NLO the 

value of A1-L from our fits is in good agreement 

MS —, '^ ^ 

with that found in Ref. g (337±28 MeV) where 
both, F2 and XF3 structure functions, have been 
fitted. There is a clear correlation between the 
effects of the NLO and NNLO approximations and 
power corrections (see table 1). At NNLO the fits 
performed with and without power corrections (in 



the IRR model) are almost equal. The significant 
decrease of the magnitude of power corrections in 
the NNLO fit with IRR model {A2 vanished within 
statistical errors) is also found with the model 
h{x)/Q^ (see Fig. 1). 

Using the value of A from the NNLO fit 
and the running of the coupling up to Af|, 
we obtain ^ as{Mz) = 0.118 ± 0.002{stat) ± 
O.QQ5{syst) ± 0.003{theory). The theoretical error 
takes into account the dependence on the initial 
Qg, the influence of the missing higher order terms 
estimated by Fade approximants and the crossing 
of the mi, threshold in the calculation of as{Mz). 

However, in the analysis there are also involved 
various approximations and shortcuts which could 
increase this uncertainty. The calculation of xF^ 
with even-n F2 anomalous dimensions P|, the 
interpolation to odd values of n |H| and the effect of 
the reconstruction method through the parameters 
a, (3 and the number of polinomials N^ax (see Eq. 
(2)) 1^,01, are not expected to affect the accuracy 
of the analysis. 

In addition, we have also studied the effect of 
using in Eq. (1) the original exponenciated formula 
for the anomalous dimension part (see Eq. (4) in 
Ref. 3). We found a change in A of 2 MeV at 
NLO and much smaller at NNLO. The effect of 
nuclear corrections has also been addressed by us 0] 
although it still deserves a more detailed study. The 
dependence with the number of active flavors (we 
work with n/ = 4) should also be carefully studied 
(see Ref. §). 

Finally, the renormalization and factorization 
scale dependence (we have fixed both equal to Q^) 
should also be estimated if one wants to make a 
meaningful precision test of pcrturbative QCD. We 
plan to present this work elsewhere pO| . 
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Table 1. Fits to Q^ > 5 GeV^ xFs CCFR'97 data. 
The starting QCD evolution point is Qq — 20 GeV^ 





Ag^ (MeV) 


^2(GeV2) 


X^ /points 


LO 


264±37 
433±53 
331±162 


-0.33±0.06 
h(x) in Fig.l 


113.1/86 

83.1/86 

66.3/86 


NLO 


339±42 
369±39 
440±183 


-0.12±0.06 
h(x) in Fig.l 


87.6/86 
82.3/86 
65.7/86 


NNLO 


326±35 
327±35 
372±133 


-0.01±0.05 
h(x) in Fig.l 


77.0/86 
76.9/86 
65.0/86 
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Figure 1. h{x) extracted from xF^ CCFR data. 
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